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Abstract: We study perturbative and non-perturbative aspects of the distribution of 
the C parameter in e + e~ annihilation using renormalon techniques. We perform an exact 
calculation of the characteristic function, corresponding to the C parameter differential 
cross section for a single off-shell gluon. We then concentrate on the two-jet region, derive 
the Borel representation of the Sudakov exponent in the large-/?o limit and compare the 
result to that of the thrust T. Analysing the exponent, we distinguish two ingredients: 
the jet function, depending on Q 2 C, summarizing the effects of collinear radiation, and a 
function describing soft emission at large angles, with momenta of order QC. The former is 
the same as for the thrust upon scaling C by 1/6, whereas the latter is different. We verify 
that the rescaled C distribution coincides with that of 1 — T to next -to-leading logarithmic 
accuracy, as predicted by Catani and Webber, and demonstrate that this relation breaks 
down beyond this order owing to soft radiation at large angles. The pattern of power 
corrections is also similar to that of the thrust: corrections appear as odd powers of A/ (QC). 
Based on the size of the renormalon ambiguity, however, the shape function is different: 
subleading power corrections for the C distribution appear to be significantly smaller than 
those for the thrust. 
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1. Introduction 

Event shape distributions have proven to be valuable in pushing forward our quantitative 
understanding of jet production and hadronization and of the interplay between perturba- 
tive and non-perturbative phenomena. Being infrared and collinear safe, these distribu- 
tions can be computed order by order in perturbation theory, without the need to introduce 
non-perturbative parameters [1-5]. Even at large center-of-mass energies (Q), however, 
event-shape distributions involve substantial power corrections. As a consequence, they 
provide an important tool to study the onset of non-perturbative physics. 

The two-jet region, where the distributions of most event shapes peak, is particularly 
important in this regard, and consequently it is also quite challenging. First of all, the 
perturbative analysis involves large Sudakov logarithms (double logarithms of the event 
shape e which vanishes in the limit of two massless jets, examples being the C parameter, 
jet masses pj, and t = 1 — T, with T the thrust). These logarithms make the perturbative 
coefficients of the distributions diverge at any finite order in the limit e — > 0. It is only upon 
resummation that one recovers the qualitative features of the cross section, namely the fact 
that it vanishes in the two-jet limit where radiation is inhibited (Sudakov suppression). A 
further complication is that power corrections, which are associated with large-angle soft 
emission, are also enhanced in this limit, as the relevant scale is typically Qe. 

The resummation of large perturbative corrections [6-15] and the parameterization 
of power corrections based on renormalon techniques [16-38] have opened up the way for 
quantitative predictions in the two-jet region. While much progress was made during the 
LEP era, some of the fundamental questions concerning power corrections have not yet been 
fully answered. A classical example is the relation ("universality") between corrections 
to different observables which may be deduced from renormalon models. Even a more 
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pragmatic motivation to study event shapes, namely to provide a precise measurement of 
the strong coupling, is hampered by theoretical uncertainties [34,36,37,39]. Thus, in spite of 
having no active e + e~ collider, further theoretical progress is important. Progress is made 
nowadays in fixed-order calculations [40], as well as in resummation and parametrization 
of power corrections [41-46]. Tools developed in the context of event shapes in e + e~ 
annihilation are often used in other applications. 

One of the aspects where the study of event shapes has taught us general lessons about 
QCD is the interplay between perturbative and non-perturbative corrections. The phe- 
nomenological success of renormalon-based models for power corrections has important 
consequences. On the one hand, it shows that perturbative tools are quite powerful. In a 
sense, the state-of-the-art theoretical description of event-shape distributions pushes per- 
turbative calculations beyond their natural regime of applicability. In general, this relies on 
the understanding that hadronization is a soft phenomenon, which does not involve signif- 
icant momentum flow and thus does not change drastically the perturbative distribution. 
On the other hand, the success of these models implies that there is no way to quan- 
tify non-perturbative corrections and relate them to matrix elements without controlling 
first perturbative corrections to all orders. This stands in sharp contrast to the standard 
approach taken when vacuum condensates are estimated in the framework of QCD sum 
rules [47]. The current understanding of power corrections highlights the significance of 
running-coupling effects at the perturbative level. 

Our main focus in this paper is the two-jet region. Quantitative predictions in this 
region are highly non-trivial, involving large logarithmic corrections from multiple soft and 
collinear gluon radiation, as well as significant running-coupling effects. The latter have 
both perturbative and non-perturbative aspects, and some of the difficulty is related to the 
ambiguous separation between them. An important feature of the two-jet region is the fact 
that, as e — > 0, it becomes necessary to resum not only singular perturbative contributions, 
but power corrections as well: in fact for values of e = 0(A/Q), which are not far from the 
peak of the distribution in typical LEP data, all power corrections of the form (A/(Qe)) p 
become equally important. The feasibility in principle of such a resummation was shown 
in Refs. [21,22], extending the factorization of the cross section in the two-jet limit to 
power-suppressed effects. Power corrections of this kind can be summarized in a "shape 
function", which can be naturally combined with the effects of Sudakov resummation. 
Renormalon calculations can then be used to construct highly constrained QCD models 
for the shape function. In this paper, we will essentially construct such a model for the C 
parameter distribution, following the Dressed Gluon Exponentiation (DGE) approach, first 
developed in connection with the thrust distribution [36]. In this approach one computes 
and resums Sudakov logs and renormalons simultaneously, obtaining information on leading 
as well as subleading power corrections, which can be summarized if desired in an ansatz 
for the shape function. Note that the combination of Sudakov effects with parametrically 
enhanced power corrections is definitely not unique to event-shape distributions. It appears 
whenever differential cross sections in QCD are evaluated near a kinematic threshold; other 
important examples are Drell-Yan production near the energy threshold [19,32], structure 
functions near the elastic limit (Bjorken x close to 1) [48,49], and fragmentation functions 
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of light [30,38] and heavy quarks [50] near z = 1. In spite of the different nature of these 
processes, certain characteristics of the perturbative expansion are generic [38], and similar 
techniques are applicable in all cases. 

As announced, we will consider here the distribution of a specific event-shape variable, 
the C parameter. This variable was first proposed 25 years ago [2], and became one of 
the classical examples of infrared and collinear safe observables. A next-to-leading order 
(NLO) 0{a 2 s ) calculation for the distribution of the C parameter was performed long 
ago [5] (the next order is still not available), while the resummation of Sudakov logs to 
next-to-leading logarithmic accuracy (NLL) has become available only recently [13]. In 
the last few years the example of the C parameter, and in particular its average value, 
had an important place in the ongoing debate concerning universality of power corrections, 
see e.g. [13,23,30,31]. On the other hand, as far as the resummation of running-coupling 
effects is concerned, the C parameter distribution was not analysed, and the prime example 
has always been the thrust [34,35]. 

The study of the thrust and of the heavy-jet mass distributions by means of DGE [36, 
37] highlighted several significant features of higher-order perturbative as well as non- 
perturbative contributions. In particular, it was found that subleading Sudakov logs, that 
are usually neglected, can give substantial contributions: they carry the factorial enhance- 
ment of renormalons. Phenomenologically, the refinement of the calculated distribution by 
DGE (as compared to the standard resummation to NLL accuracy) turned out to have a 
dramatic effect on data fits in the peak region. The most impressive demonstration of this 
fact is the agreement between the non-perturbative parameters extracted from the thrust 
and the heavy-jet mass distributions [37]. The effect this resummation has on the ex- 
tracted value of the coupling is also substantial. It is, of course, of great interest to extend 
the DGE analysis to other event-shape variables and learn which features are generic and 
which depend on the observable. The purpose of this paper is to perform such an analysis 
for the C parameter. 

We proceed as follows. In the next section we recall the definition of the C parameter, 
present the kinematics in a three-particle final state where the gluon is off shell, and 
compute the characteristic function. The final result is summarized in Appendix A. In 
Section 3 we expand the characteristic function at small C (details are given in Appendix 
B) in order to identify the source of logarithmically enhanced terms which dominate the 
two-jet limit; next, we construct a Borel representation of the Sudakov exponent, and we 
use it to extract perturbative as well as non-perturbative information on the distribution, 
comparing it with the case of the thrust. Our results are summarized and briefly discussed 
in Section 4. 

2. The characteristic function 

The C parameter for electron-positron annihilation events was originally 1 defined [3, 4] as 

C = 3(A 1 A 2 + A 2 A3 + A 3 Ai) , (2.1) 

1 Another definition has been introduced in [2]; see also [5]. 
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where X a are the eigenvalues of the matrix 




1 




(2.2) 



and Pa are the spatial components (a = 1,2,3) of the i— th particle momentum in the 
center-of-mass frame. The sums over i and j run over all the final state particles. 
A related definition in terms of Lorentz invariants is 



a discussion of mass effects in power corrections, see [51]). 

The C parameter varies in the range < C < 1. C = corresponds to a perfect 
two-jet event (with massless jets), while C = 1 characterizes a spherical event. Planar 
events, including, in particular, the 0(a s ) perturbative result, are distributed in the range 
< C < 3/4. 

The information needed to perform renormalon resummation at the level of a single 
dressed gluon (the large-/3o limit) and, eventually, learn about power corrections, is present 
in the leading-order differential cross section, calculated with an off-shell gluon [16,27,52], 
sometimes called "characteristic function" . In this section we present an exact calculation 
of the characteristic function for the distribution of the C parameter. An analogous cal- 
culation for the thrust was performed in [35]. Although the calculation and the general 
structure of the result are similar, in the case under consideration the expressions involved 
are significantly more complicated. 

It should be noted that the renormalon calculation we perform treats the decay of the 
gluon inclusively. The characteristic function depends only on the total invariant mass 
of the particles eventually produced by the gluon and not on their separate momenta. 
Since event-shape variables such as the C parameter are sensitive to the momenta of 
individual final-state particles, our calculation differs, for example, from a strict large-iVy 
limit. This point was first noted in [25], and was since addressed in various occasions, 
first in applications to single-particle inclusive cross sections [33], then in the context of 
the conjectured universality of power corrections [29]. In the thrust case the inclusive 
approximation is good, in the sense that higher-order perturbative terms are numerically 
close to the strict large-A^ result [34-36] (larger differences occur if one considers purely 
non-Abelian contributions). The approximation is especially good in the two-jet limit in 
which we are interested. In the case of the average C parameter, Smye [31] has performed 
a detailed analysis of the effect non-inclusive contributions have on the coefficient of the 
1/Q correction, within the framework of [29]. Also in this case the inclusive approximation 
is close to the large-Nf result, while larger differences appear in the non-Abelian part. An 
analysis of this kind has not yet been performed for the distribution. We observe also that 
it is unclear to what extent non-inclusive contributions can consistently be treated within 
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where is the four-momentum of the i-th particle and q is the total four momentum, 
q = ■ The two definitions are equivalent provided particle masses are neglected (for 
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the framework of the dispersive approach, since they contain terms that are genuinely 
unrelated to running-coupling effects, even in the abelian limit. 

It should be kept in mind that, with present technology, renormalon models can only 
be used as a tool to gather information on the general structure of power corrections, and 
possibly a rough estimate on their size. For such purposes, we believe that the inclusive 
approximation is sufficient, a view which is supported by the phenomenological success of 
the data fits performed with it, and also by the small relative size of the non-inclusive 
correction in the cases in which it was computed. Given the complexity of the calculation 
in the case of the C parameter, we will not attempt here a calculation in the strict large- Nf 
limit. A calculation of the distribution along the lines of [31] would in any case be welcome, 
since it would improve the control on phase-space effects, and gauge the stability of our 
result. 

In order to perform the calculation of the characteristic function we first calculate the 
C parameter for a three-particle final state: a quark and an antiquark with momenta p\ 
and P2 (pi = p\ = 0) and an off-shell gluon with momentum k. It is convenient to define 
the normalized gluon virtuality by £ = k 2 /q 2 , and the variables 

xi = 2pi • q/q 2 , 

x 2 = 2p 2 ■ q/q 2 , (2.4) 
y = 2k ■ q/q 2 = 2 - xi - x 2 , 

which correspond to the energy fractions in the centre-of-mass frame. Using these variables 
one finds 

I n - nit (l-x 1 )(l-x 2 )(l-j/ + 2g)-£ 2 

c(xi,x 2 ,£) = C 6 = . (2.5) 

xix 2 y 

Here we rescaled the shape variable by a factor of 6; we will compute the differential 
cross section for the rescaled variable c, from which the standard observable can be readily 
obtained, as 

-^(C,Q 2 )= 1 - 1 -^(c,Q 2 ) . (2.6) 

a aL o a dc c =c/6 

Depending on the precise interpretation of the definition of the C parameter for massive 
particles there can be in Eq. (2.5) an additional term (from i and j in (2.3) both corre- 
sponding to the gluon) of the form — £, 2 /y 2 . Here, however, we will be using Eq. (2.5) as 
written. As we shall see below, in the small-c region, c = 0(£/y), so this term is negligible. 

The renormalon-resummed differential cross section in the single dressed gluon ap- 
proximation is 

Ida 2 \ Cf f 1 d.F(£,c) Aler .2\ 



,-(c,Q i ) = -^f / d^^^A^Q 2 ), (2.7) 
cr dc 2fj J d£ 

where /3q = — ^Nf, and A(£Q 2 ) is the large-/?o running coupling (A = (3q0. s /it) on 

the time-like axis, admitting the Borel representation 

A(iQ 2 ) = Hdu (Q 2 /A 2 y u ™et«r M , (2.8) 
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where A is the QCD scale in the MS scheme, the constant 5/3 comes from the renormal- 
ization of the fermion loop in this scheme, and the sine factor originates in the analytic 
continuation to the time-like axis. 

The characteristic function is of the form 



dxidx 2 M(xi,X2,ti)5 (c(xi,x 2 ,0 - c) 



(2.9) 



Here M is the squared matrix element for 7* — s 
factor extracted), and is given by 

(l-xi)(l-x 2 ) (1-xi) 2 {I-X2) 2 
Phase space (illustrated in figure 1) is limited by the restrictions 



qqg (with the coupling and the colour 



(2.10) 



hard limit 
soft limit 



X\ + x 2 > 1 - £ , 
{l- Xl )(l-x 2 )>S 



(2.11) 



where the hard and soft limits correspond to the upper and lower bounds on the gluon 
energy fraction y = 2 — x\ — x 2 . 

Note that along the soft limit, which will be our main interest in the next section, the 
expression for c in Eq. (2.5) simplifies a lot. One gets 



c{xi,x 2 ,o\soft = Z/y ■ 



(2.12) 



Along this phase-space boundary, the maximal value for c with fixed £ is obtained for 
x\ = x 2 = 1 — y/£, so y = y m i n = 2\/£ and c^x = This is the relevant phase- 

space limit for large-angle soft gluon emission. The minimal value of c is obtained in the 
intersection with the hard limit, for either x\ = or x 2 = 0, so y = y max = 1 + £ and 
Cmin = +0- This is the relevant limit for collinear gluon emission. The maximum 
value of c is attained away from the boundaries of phase space, as seen in Figure 1. It 
can be determined by maximizing the explicit expression in Eq. (2.5) within the physical 
region, and for small £ is given by c max = 1/8 + 3£/4 + 0(£ 2 ). 

In order to perform the integrations in (2.9), it is convenient to change variables from 
x\ and x 2 into y and z = (1 — x 2 )/y. In the new variables the integral has the form 



r^r^^KzsUz-z^z-z-) 

J2V1 y Ji(i-V^w) V 



y{l + 2i-y{l + c)) 

( i_ y (i_ z ))2 + (1 _ yz) 2 + ^ 2 +i-y) e e 



z(l-z) 



z 2 {l-z) 2 \ ' 



where 



z± 



1 ± 4/I 



A(cy(l-y) + e) 
y 2 (l + 2e-y(l + c)) 



(2.13) 



(2.14) 



Here the symmetry x\ 



x 2 appears as z 



1 — z. Using this symmetry, the integral 



over z equals twice the integral between the lower limit in Eq. (2.13) and z = 1/2, where 
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Figure 1: The x\ - X2 plane for a fixed £ = 0.06. The phase-space limits (2.11) (green) are shown 
together with contours of fixed c belonging to the two regions in (2.18): contours that close within 
the phase-space boundaries, corresponding to Fh (thick, blue), and contours that reach the soft 
phase-space limit, corresponding to Fi (thin, red). 



only z = Z- in the 5 function is relevant. The condition z_ > | ^1 — y/l — 4£/y 2 j for 
physical values of c implies that 

- < y < = y . (2.15) 

c~l+c v/ 

The condition that z_ be real, on the other hand, implies that 

(1 + c)y 3 - (1 + 2£ + 4c)y 2 + 4cy + 4^ 2 < . (2.16) 

It turns out that within the physical region all three roots 2/1,2,3 of Eq. (2.16) are real. 
Choosing y% < 7/2 < 2/3, Eq. (2.16) requires that 

y 2 < y < y 3 ■ (2.17) 
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Investigating the solutions of Eq. (2.16), one finds that 1/2 > £/c for c > \/£/2, while 
V2 < £/c f° r c < "v/f/2, so the conditions (2.17) and (2.15) together imply that the lower 
integration limit over y is 1/2 for c > \/f/2, but it is £/c for c < y/%/2. The result is 
therefore 



^i(C.c) = / dyf(y,c,0 

h/c 

rvz 

Fh(£,,c) = / dyf(y,c,0 

J V2 



i+§ — — 2 ' 



(2.18) 



2 



< c < c max = I + f£ + 0(£ 2 ) , 



where 



/(y,c,0 



y(l + ^-y) s 



(1 + c) 2 ( cy (i - y ) + e 2 ) 2 y/(y -y) 3 {y 3 -y)(y-yi)(y-y2) 

x - c (1 + c) y 4 + ((5 c + 2 c 2 + 1)£ + 3 c + 4 c 2 ) y 3 

+ ((-2 -4c) £ 2 + (-2c 2 - 1 - 10cK -4c-4c 2 )y 2 
+ ((-2c-2K 3 -2£ 2 +8£c + 2c)2/ + 4£ 4 + 8£ 3 + 2£ 21 



(2.19) 



The two regions of the parameter space of Eq. (2.18) are distinguished by whether, for 
fixed £ and c, the soft phase-space limit in Eq. (2.11) can be reached. The corresponding 
classes of contours are shown in Figure 1 in the X\ - X2 plane. A similar situation occurs 
in the case of the thrust distribution [35] . In both cases the characteristic function is given 
by two different analytic function in these two regions. 

Finally, c) can be expressed as a sum of standard elliptic integrals. These are com- 
plete elliptic integrals for Fh{€i c ) an d incomplete ones for JF^(£, c). The explicit expressions 
are given in Appendix A. 

3. The Sudakov exponent 

Having calculated the characteristic function for the distribution of the C parameter, we can 
use DGE to compute the Sudakov exponent for the two-jet limit in the large-/?o limit. The 
procedure we follow was introduced in [36], where the exponent for the thrust distribution 
was computed. Similar calculations were done since for heavy-jet mass distribution [37], 
light-quark fragmentation, deep inelastic structure functions and Drell-Yan production for 
x — > 1 [38], and, most recently, the heavy-quark fragmentation function [50]. 

The Borel representation of Eq. (2.7) can be constructed by using the explicit expres- 
sion for A(£Q 2 ) given in Eq. (2.8) and changing the order of integration. One obtains 



1 da 

a dc 



(c,Q 2 ) 



SDG 



= W / du(Q>/Ar U B(c,u), 
Z P0 Jo 



where the subscript stands for a Single Dressed Gluon. The Borel function is 



B(c, u) 



smiru 5, 
es' 



TTU 



'lc 2 



dFh{i,c) 



d£ 



rc/(i-c) ^ dmjC) 

/4c 2 



J 4< 



da 



(3.1) 



(3.2) 
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While these expressions are completely general, and can be used to resum the renor- 
malons in the large—/3o limit for arbitrary c, we are interested specifically in Sudakov logs, 
which dominate the cross section in the two-jet region. Starting from the full character- 
istic function, the first step is to identify the 0(l/c) singular terms, namely those terms 
that upon integration over the gluon virtuality in Eq. (2.7) or in Eq. (3.2) would lead to 
logarithmically enhanced contributions in the perturbative expansion of the distribution. 

The identification of the relevant terms was straightforward in the case of the thrust, 
but it is much less so here, given the complexity of the characteristic function, Eq. (2.18). 
One can expand the integrand given in Eq. (2.19) or the closed-form expressions of Ap- 
pendix A at small c, but this requires of course to specify how £ behaves in the limit 
considered. Given the limits of integration in Eq. (3.2), there are two natural c — > limits 
to consider: one with fixed a = £/c, which is relevant in the collinear region, and the second 
with fixed b = £/c 2 , relevant in the large-angle soft emission region. 

On general grounds, one expects Sudakov logs to emerge from the soft boundary of 
phase space in Eq. (2.11). Figure 1 shows that only for T\ the soft boundary of phase 
space is relevant, so it is only the second term in Eq. (3.2) that is expected to be relevant. 
We will verify this statement explicitly below. A similar situation occurs for the thrust 
distribution [36]. 

The details of the small-c expansion of J-\ at fixed a and at fixed b, and of at fixed 
b, are summarized in Appendix B. The final results for the leading 0(l/c) terms are 

^(e,c)| fixcda ^ /c = \ [-4 In a - 3 + 2a + a 2 ] + O(c ) 

•^(£> c )lfixcd6=g/c2 = ~ 

while 

MZ,c)\ &xcdb ^ /c2 = - [-41nc - 3] + O(c ) . (3.4) 

Eq. (3.4)) implies that indeed the first term in Eq. (3.2), the one involving Th, is irrelevant 
for logarithmically enhanced terms. In order to compute the Sudakov exponent we therefore 
concentrate on F\. Using Eq. (3.3) it is straightforward to write down an expression for F\ 
that has the correct 0(1/ c) terms both for fixed a and for fixed b. It is given by 



41n6c-3-8hW - (l + y/l-4/b) 



+ 0(c°), (3.3) 



•Fj(£,c)| 0(1/c) 



- 1 in I - 



Taking a derivative of Eq. (3.5) and substituting into Eq. (3.2) we find 



. (3.5) 



B(c, u) 



5,, smiru 
2es w 



TTU J4 C 2 



J Ac 



l £ 



+ 



8c 



(3.6) 

where we simplified the upper limit without affecting the relevant (log-enhanced) terms. 
As expected from Eq. (3.3), only the first and the last terms in the square brackets in 
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Eq. (3.6) contribute to the Sudakov exponent owing to the lower (large-angle) integration 
limit, while the first three terms, and not the last, contribute owing to the upper (collinear) 
limit. The result of the integration is 



B(c, u) 



logs 



2e3 u 

iru 



4(2c) 



-l-2u 



— c 



-l-u 



2 

u + T 



u 



+ 



1 



• (3-7) 



This simple expression summarizes all log-enhanced contributions to the c distribution 
(3.1), to any order in perturbation theory, in the large-/?o limit. The result is written as 
a sum of two ingredients, each originating in a different region of phase space: the first, 
where the relevant scale is AQ 2 c 2 (gluon virtuality £ ~ 4c 2 ) is associated with large-angle 
soft emission (it can be traced through the calculation to the region where y ~ 2y / £ and 
x\ ~ X2 ~ 1 — \/£); tli e second, where the scale is Q 2 c (gluon virtuality £ ~ c) is associated 
with collinear emission (corresponding to y ~ 1 + £, and either x\ ~ or 22 ^ 0). Note 
that infrared safety is guaranteed in Eq. (3.7) thanks to the cancellation between 0(l/u) 
terms between the soft and the collinear ingredients, so each of these ingredients is ill- 
defined unless a factorization prescription is introduced. Note also that Eq. (3.7) is free of 
renormalon singularities, but the Borel integral (3.1) does have convergence constraints at 
small c, which will turn into renormalons in the Sudakov exponent upon performing the 
necessary Laplace transform, as outlined below. 

For comparison we quote the analogous result for the thrust distribution [36,38]. Defin- 
ing t = 1 — T, the Borel function in the large-/?o limit is 



B(t,u)\ logs = 2es« 



sin TTU 



TTU 



j.—l—2u _ j—\—u 



U 



2 1 1 
-+ 1 + ~ 

u 1 — u 2 — u 



(3.8) 



The result for the heavy-jet mass [37] (in this approximation) is the same as Eq. (3.8) up 
to an overall factor of 1/2. The similarity, as well as the differences, between Eq. (3.7) and 
Eq. (3.8) are easily understood: the collinear ingredient is identical (with c replaced by t), 
while the large-angle soft emission ingredient is different. In fact, replacing c by 1 — x, the 
collinear ingredient coincides with the large-/3o expression [38] for the jet function [7-9] 
which controls the large-x limit in single-particle inclusive cross-sections in e + e~ annihi- 
lation [38,50], as well as for structure functions near the elastic limit [48,49]. This object 
describes the radiation from an unobserved jet under a restriction on its invariant mass. It 
therefore appears in each of these observables. 

As noted above, the large-angle soft emission ingredient for the c distribution in 
Eq. (3.7) turns out to be quite different from that of the thrust, in Eq. (3.8). This is 
not surprising: different event-shape variables weigh differently the soft radiation. One 
immediately recognizes that the relevant scale of large-angle soft radiation in the c distri- 
bution is 2Qc, while for thrust it is Qt. As we shall see below this factor of 2 suppresses 
subleading logarithmic corrections, as well as subleading power corrections, in the case of 
c, as compared to t. 

We now turn to the computation of the Sudakov exponent. Owing to the additivity 
property of the C parameter in the small C limit with respect to multiple emission, and the 
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factorization property of QCD matrix elements for soft and collinear radiation, logarithmi- 
cally enhanced terms in the perturbative expansion exponentiate in Laplace space. Conse- 
quently, the resummed cross section can be written, in full analogy with the thrust [10,36], 
as 

rk+ioo j 

= / e » c exp[S(v,Q 2 )] , (3.9) 

DGE Jfe-ioo Z7Tl 

where k is to the right of the singularities of the integrand. The Sudakov exponent is 
given by 



a ac 



S(u,Q 2 )^(e-^)= r dc^(c,Q 2 ] 
Jo ° dc 



(e-- - 1) , (3.10) 



where the upper limit of integration was extended to infinity (i.e. beyond the range where 
the physical distribution has support) to comply with the standard definition of the Laplace 
transform. The integral is dominated by the Sudakov region, so that singular terms for 
c — > in Eq. (3.1) generate logarithmically divergent terms for v — > oo in S(v,Q 2 ). It 
should be noted that, as far as these singular terms are concerned, and to any logarithmic 
accuracy, the Laplace transform is equivalent to a Mellin transform, where v is the moment 
conjugate to 1 — c. 

Next, we compute S(u,Q 2 ) in the large-/3o limit using Eq. (3.1) and Eq. (3.7). The 
result is 

poo 



S(u,Q 2 ) = / du (Q 2 /A 2 P Bl(u) , (3.11) 
l Po Jo 



with 



Bl(u)= r d cB(c,u)(e-^-l) = 2e^ S ^^ T(-2u) (v 2u - l) 2 l ~ 2u 

Jo v ; 7TU [ ; T(i + u) 

- r(-u) {v u -\)(- + — L_ + J-) ] , (3.12) 
\u 1 — u 2 — u J 

where terms suppressed by powers of \/v were neglected. 

Contrary to Eq. (3.7), the Borel function in Laplace space given by Eq. (3.12) has 
renormalon singularities. The appearance of infrared renormalons as a consequence of 
taking a Laplace or a Mellin transform is characteristic of differential cross-sections near a 
kinematic threshold [36,38,49,50]. In Eq. (3.12) renormalons appear in both the large-angle 
soft emission factor and in the jet function. In the former, they correspond to corrections 
that scale as odd powers of Q/is, while in the latter to ones that scale as the first two 
powers of Q 2 /is, i.e. twist four and twist six, as in structure functions. The leading power 
corrections in the Sudakov region (so long as Q 2 c 3> A 2 ) are of the former type, and they 
can be summed up into a shape function, as discussed in [19-23,36,37]. We further address 
this issue below. 

Eq. (3.12), together with Eq. (3.11) and Eq. (3.9), is our final result for the c dis- 
tribution, calculated by DGE. The experience gained in the case of the thrust and the 
heavy-jet mass distributions [36, 37] shows that there is a significant difference between 
the resummed cross section a la DGE and standard resummation with NLL accuracy. Ow- 
ing to the renormalon factorial growth, the additional terms (subleading logs) which are 
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resummed by DGE and neglected otherwise are numerically important. DGE has a sig- 
nificant impact on the quality of the description of the distribution in the peak region, on 
the extracted non-perturbative parameters, and, in particular, on the consistency of the 
latter between the thrust and the heavy-jet mass distributions. It is therefore important 
to analyse more event-shape distributions with the same methodology, and the road is now 
open to do so for the C parameter. 

It should be emphasized that controlling the large-/?o terms is sufficient for generating 
the entire set of logs up to NLL accuracy. This can be done [36] by promoting the running 
coupling of Eq. (2.8) to two loops, while replacing the constant 5/3 in the exponent by the 
full coefficient of the x — ► 1 singular term in the NLO splitting function [12], according to 

One should keep in mind, though, that this replacement is sufficient only to NLL accuracy, 
while the correct generalization of the Borel function (3.12) beyond the large-/3o limit is 
not known. See [49] for a relevant discussion. 

Of course, as it stands, the integral in Eq. (3.11) is ill defined due to its renormalons 
singularities. It is straightforward, however, to perform a principal-value regularization 
either analytically [36,37] or numerically [50]. The size of the ambiguity of Eq. (3.11), as 
measured by the residues of the singularities, can be used as an estimate of the size of the 
corresponding power corrections, as discussed below. Let us also recall that, in order to use 
the DGE result in practice, one needs to match it onto the known fixed-order calculation, 
to take into account terms which are not logarithmically enhanced. The procedure [36] is 
similar to the one used in a standard NLL resummation program [10]. 

The location of the renormalons in Eq. (3.12) is identical to the case of the thrust. For 
easy comparison we quote the analogous expression derived in that case, where v is the 
Laplace conjugate variable to t = 1 — T. It is given by 

r(_2«) - 1) - - r(-u) K + + J-) 

' u \u 1 — u 2 — u J 

(3.14) 

In order to extract the perturbative coefficients of the exponent S(u, Q 2 ) in the large-/3o 
limit, it is sufficient to expand B v {u) in powers of u and replace u n by n\ (Poa s /Tr) n+1 . 
Before looking at the actual coefficients, let us first note that the difference between the 
c distribution, Eq. (3.12), and the thrust distribution, Eq. (3.14), which is all due to the 
coefficient of {y 2u — l) , is relatively small. Taking the ratio between these two coefficients 
and expanding in u we get 

2 -^M = 1 -^ + 2C3 „3 + 0(l( 4 ) . (8JS) 

r(2+u) 6 

The fact that the ratio is 1 at leading order (as it must be, in order to cancel the l/u 
singularity of the jet function, which is identical in the two cases) implies that the leading 
logs (L n+1 a™ where L = Inv, for any n) are identical for the two shape variables; the 
fact that the next term, of order u, is missing, implies that also the next-to-leading logs 



Tit/ \ o sm7ru 
Bl(u) = 2e3 u 



TTU 
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(L n a™, for any n) are identical. This confirms the prediction by Catani and Webber [13]. 
Differences between the two distributions appear only at the next-to-next-to-leading log- 
arithmic level. 

In order to illustrate the enhancement of subleading logarithms [36], as well as the 
difference between the c distribution and the thrust in this context, we present below an 
expansion of Eq. (3.12) and Eq. (3.14) to the first few orders. One finds for the c parameter 

B c u (u) = —21? + 0.691 L 

+ (-2 L 3 — 5.297 L 2 + 0.095 L) u (3.16) 

+ (-1.167 L 4 - 5.527 L 3 - 7.911 L 2 - 22.709 L) u 2 

+ (-0.5 L 5 - 3.262 L 4 - 7.943 L 3 - 30.058 L 2 - 61.322 L) v 3 

+ (-0.172 L 6 - 1.405 L 5 - 4.639 L 4 - 22.488 L 3 - 67.59 L 2 - 135.22 L) u 4 

while the thrust gives 

^thrust^ = _2 L 2 + 0.691 L 

+ (-2L 3 - 5.297 L 2 - 6.485 L) u (3.17) 

+ (-1. 167 L 4 - 5.527 L 3 - 14.491 L 2 - 31.655 L) u 2 

+ (-0.5 L 5 - 3.262 L 4 - 12.329 L 3 - 39.003 L 2 - 80.940 L) u 3 

+ (-0.172 L 6 - 1.405 L 5 - 6.832 L 4 - 28.452 L 3 - 87.21 L 2 - 175.80 L) n 4 

As explained above, the leading and next -to-leading logarithms are the same. Beyond this 
order the general trend is similar, however the increase of the coefficients is milder in the 
case of c. 

Finally, we return to the issue of power corrections. One can construct a parametriza- 
tion of non-perturbative corrections based on the ambiguity induced by renormalon singu- 
larities in B v {u). The residue of a pole at u = m/2, where m is an odd integer, multiplies 
an ambiguous contribution, and thus a non-perturbative correction of order (Ai//Q) m in 
S(u, Q 2 ) is necessary to compensate for the ambiguity of the perturbative result. Summing 
over m, these corrections amount to a multiplication of the perturbative Laplace-space re- 
sult exp [S(u, Q 2 )] , entering Eq. (3.9), by a non-perturbative shape function [19-23,36,37] 
of the variable Au/Q. Since these corrections exponentiate, the leading m = 1 correction 
generates a shift of the perturbative distribution [20,26], while m = 3 and higher-order 
powers generate smearing. Although the large-/?o renormalon calculation is not sufficient 
to determine the actual magnitude of power corrections, the size of the residues 2 may be 
taken as a naive estimate. The first few residues (times tt) are summarized in Table 1. 

Assuming that the large-/?o residues do provide some hint on the size of the corrections, 
from Table 1 one would conclude that the non-perturbative shape functions for the c 
parameter and the thrust distributions must be significantly different. From the perspective 
of this renormalon model, it seems therefore that the conjecture of [23], that the same 

2 Note that the large-/3o limit does not necessarily provide a good estimate of the residue. Moreover, in 
the full theory the singularities are usually not simple poles. 
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Correction 


residue (c) 


residue (t) 


ratio (c/t) 


{Kv/qY 


2vr 


8 


f ~ 0.79 


(Xv/Qf 


_ 72 


4 
27 


|f ~ 0.29 


(Au/Q) 5 


7T 

12800 


1 

375 


Mf - °- 092 


(Au/Q) 7 


7T 

3612672 


1 

30870 


4096 — U - UZ ' 


(Au/Q) 9 




1 


315 71- ~ n nnvg 


1528823808 


3674160 


131072 — u.uuiu 



Table 1: The size of the residues of rcnormalon singularities in the large-/?o limit, based on 
Eq. (3.12) for the c distribution and on Eq. (3.14) for the thrust t. The numbers quoted are ir 
times the coefficients of (Au/Q) , where we define A 2 = A 2 e 5 / 3 . We ignore here the 0(1) factor 
C F /(2p ) in Eq. (3.11). 

function (with the same parameters) would be appropriate for both distribution, should 
be excluded. In fact, while the powers of Au/Q identifying the relevant moments of this 
function are the same for both variables, subleading corrections to the c distribution are 
significantly smaller than the corresponding ones for the thrust. This means that while the 
shift of the two distributions is of comparable magnitude (the first line in the table) , the 
extent to which subleading non-perturbative corrections smear the c distribution is much 
smaller than for thrust. 

4. Conclusions 

In this paper we performed a renormalon calculation for the distribution of the C param- 
eter. We adopted the dispersive approach and extracted all-order information from the 
differential cross section with a single off-shell gluon. The only approximation we made 
with respect to the exact, all-order result in the large-iV/ limit, was to treat the gluon 
decay inclusively. 

The characteristic function given in Appendix A can be used to improve existing fixed- 
order calculations of the distribution and its first few moments through the resummation 
of running-coupling effects [34,35]. We recall that the perturbative coefficients of the first 
few moments of event-shape distributions, and in particular the average, are dominated 
by running-coupling effects, so the impact this resummation has on phenomenology is 
significant. 

Here we concentrated on the distribution in the two-jet region, where the perturbative 
expansion is dominated by Sudakov logs. Starting from the exact characteristic function, we 
identified the origin of logarithmically enhanced terms and computed the Sudakov exponent 
in the large-/?o limit by Dressed Gluon Exponentiation, similarly to previous calculations 
for the thrust and the heavy-jet mass distributions [36,37]. We showed that the all-order 
result for the Sudakov exponent, given in Eq. (3.12), separates into two ingredients in a 
natural way: one is the jet function, depending on Q 2 c, where c = C/6, and the other 
is associated with soft emission at large angles, with momenta of order 2Qc. The jet 
function [7-9] describes the radiation from an unobserved jet under a restriction on its 
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invariant mass. It can be denned and computed in a process independent way, and it plays 
a role in a large class of differential cross sections near partonic threshold, which include, 
in addition to event-shape distributions [36,37], the coefficient functions for single-particle 
inclusive cross sections in e + e~ annihilation, as well as deep-inelastic structure functions. 
The same large-/?o result for the jet function was shown to be relevant for all of these 
observables [38,48-50]. 

Throughout the paper we compared the case of the C parameter to that of the thrust. 
As noted first by Catani and Webber [13], the two distributions are closely related: upon 
scaling C by 1/6 the Sudakov exponents of the two variables coincide to NLL accuracy. 
Since the jet functions are the same to any logarithmic accuracy, differences between the 
two exponents appear only in the large-angle soft emission ingredient. The large-/?o results, 
Eq. (3.12) and Eq. (3.14), indicate that such differences appear at the NNLL order. The 
universal phenomenon that subleading logs appear with increasing numerical coefficients 
owing to infrared renormalons [36] is realised, of course, also for the C parameter. The 
increase of the coefficients, however, is milder than for thrust. The difference is largest for 
the NNLL term at order 3 0(u), i.e. 0(a 2 s ), in Eq. (3.16), as compared to Eq. (3.17). 

It is natural to expect that the difference between the thrust and the c distribution, 
owing to large-angle soft emission, would be realised also at the non-perturbative level. 
While the general pattern of renormalon singularities in the Sudakov exponent is similar in 
the two cases, and the parametrization of non-perturbative corrections as a shift [20,26] of 
the perturbative distribution, or better, through a convolution with a shape function [19- 
23,36,37], is appropriate in both, it seems that the corrections themselves are different. 
This conclusion is based on the comparison (Table 1) between the renormalon residues 
in the two cases: subleading power corrections for the c distribution are expected to be 
smaller and thus a shift should be a better approximation for the c distribution than it is for 
thrust. We emphasize that the residues computed in the large-/3o limit are not necessarily 
indicative of the actual size of the corrections, and it remains to be seen whether future 
phenomenological analyses will support these conclusions. 
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3 Note that the coefficients at this order coincide with the exact large Nf result - they are not influenced 
at all by the inclusive approximation [36,37]. 
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A. Explicit expression of the characteristic function 



The last integration over the gluon energy fraction y in Eq. (2.18) can be explicitly per- 
formed, and the result expressed in terms of standard elliptic functions. This is a con- 
sequence of a well-known theorem of Legendre, stating that any integral of the general 
form 

l(x) = j dxR x,S(x) , (A.l) 

where R is a rational function of its arguments, while S 2 (x) is a polynomial of degree d < 4, 
can be expressed as a linear combination of the three basic elliptic integrals, plus the integral 
of a rational function. Our integrand, Eq. (2.19), clearly fulfills the requirements. To fix 
our notation, we define the three basic kinds of (incomplete) elliptic integrals by 



(0,m) ee / d6 = / 

J o \/l - ms'm z 9 Jo 



dx 



\J\-m sin 2 9 Jo a/(1 - x 2 )(l - mx 



2\ 



r<f> /-sin^ R _ 2 

E{<f>,m)= / deVl-msm 2 9 = / J , (A.2) 

Jo Jo V 1 — x 

r sin0 dx 



U(n,<p,m) ee / d8 



Jo 



{l-nsm 2 9)Vl-msm 2 e Jo (1 - nx 2 )y/(l - x 2 )(l - mx 2 ) 



The corresponding complete elliptic integrals are obtained by setting <\> = n/2, according 
to K(m) ee F(ir/2,m), E(m) ee E(ir/2,m), and Tl(n,m) = U(n,ir/2,m). The algorithms 
to reduce the generic integral in Eq. (A.l) to a combination of the standard ones given 
in Eq. (A.2), as well as the analytic properties of elliptic integrals, are described in some 
detail in Ref. [53]. 

To express the result of the integration over the gluon energy fraction y in Eq. (2.18), 
it is convenient to recall the special values of y that correspond to singularities in the 
integrand, Eq. (2.19). First, there are the three roots of the cubic equation (2.16), which 
are real in the physical region and which we ordered according to y\ < y<i < y%. Recall that 
the y integration specified in Eq. (2.18) extends over the range from y2 to 2/3 for J-^, or part 
of it, from a = £/c to 2/3, for T\. Thus these square-root singularities are either outside 
the integration region or on its boundaries. An additional singularity in the integrand 
appears at y = yo = (1 + 2£)/(l + c), however it obeys yo > y% so it is always outside 
the integration region. Finally there are two double poles outside the integration region at 
y = y 4 = (1 - Vl + 4£ 2 /c)/2 and at y = y 5 = (1 + y/l + 4£ 2 /c)/2. 

The arguments of the relevant elliptic integrals can be expressed in terms of simple 
ratios of these special values of the gluon energy fraction. One needs 



X (£, c) = arcsin 



{yo - 2/2X3/3 - a) 

(2/3 - 2/2) (2/0 - a) 



(A.3) 



m „ c) = (2/3-2/2)(2/0-2/i) (A 4) 

(2/0-2/2X2/3-2/1) ' 

n(e,c) = ^^, (A.5) 
2/0 - 2/2 
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n 4 (£, c) 
"-5 (£, c) 



(1/3 - 2/2X2/4 - 2/0) 
(2/0 - 2/2) (2/4 - 2/3) ' 

(2/3 - 3/2X2/5 - 2/0) 



(A.6) 



(A.7) 



(2/0-2/2X2/5-2/3) ' 

In terms of these rather intricate functions of £ and c, one can write the explicit expression 
for J^(£, c). It is 

^ (£, c) = r (£, c) + / (£, c) F [ x (£, c) , m (£, c) ] + e (£, c) £ [x (£, c) , m (£, c) 

+ p (£, c) n n (£, c) , % (£, c) , m (£, c) + p 4 (£, c) n n 4 (£, c) , x (£, c) , m (£, c) 



+ P5 (£, c) n n 5 (£, c) , x (£, c) , m (£, c) 



(A.8) 



where 



r(£,c) 



4£ 3 + c£(l + 18e + 32£ 2 ) 



c 3 (i + c) 2 ( c +4e 2 ; 

+ c 2 (4 + 7£ + 20£ 2 + 48£ 3 ) + c 3 (4 + 10£ + 2£ 2 + 4£ 3 ) 1 , (A.9) 
and the coefficients of the five elliptic integrals involved are given by 

^' c ) = ^7 ^ 1 77 

c z (l + c) (c + 4£ 2 )yi X /(2/3-2/i)(2/o-y2j 

4 (1 + c) £ 2 (l2£ 3 + 2c £ (1 + 6£ + 18£ 2 ) + c 2 (3 + 6£ + 24£ 2 + 68£ 3 ) 

+ 2c 3 (3 + 7£ + 2£ 3 ) ) + 4c yi (l0£ 3 (1 + 2£) + c£ (2 + 21£ + 88£ 2 + 120£ 3 + 16£ 4 ) 
+ 2c 2 (2 + 10£ + 25£ 2 + 52£ 3 + 58£ 4 ) + c 3 (5 + 22£ + 45£ 2 + 64£ 3 + 64£ 4 + 16£ 5 ) 
+ 2c 4 (1 + if (2 + 3£) ) - y 2 (1 + c) (1 + 2£) (l2£ 3 + 2c£ (l + 6£ + 18£ 2 ) 



+ C 



(3 + 6£ + 24£ 2 + 68£ 3 ) + 2c 3 (3 + 7£ + 2£ 3 ) ) 



(A.10) 



e(£,c) = - 



\/(2/3 - 2/i) (2/0 - 2/2) 



12£ 3 + 2c£(l + 6£ + 18£ 2 ) 



c 2 (l + C ) 2 (c + 4£ 2 ) 
+ c 2 (3 + 6£ + 24£ 2 + 68£ 3 ) + 2c 3 (3 + 7£ + 2£ 3 ) 



(A.ll) 



p(£,c) 



4 (yo - 2/3) 



c 2 (1 + c) 3 ^(2/3 - 2/i) (2/0 - 2/2) 
+ c 2 (1 + 6£ + 13£ 2 ) + c 3 (2 + 6£ + 9£ 2 ) 



3£ 2 + c (2 + 4£ + 11£ 2 ) 



(A.12) 



Pa (£, c) 



32 (y - 2/3) e 2 



(2/5 - 2/4) (1 - 22/3 + 2/4 - 2/5) (1 - c + (1 + c) (y 5 - 2/4) + 4£) 3 (c + 4£ 2 



-17- 



12£ 7 + 2c£ 5 (5 + 34£ + 172£ 2 ) 



c 5 v(y3 - vi) (yo - V2) 

+ c 2 £ 3 (1 + 80 (2 + 37£ + 18£ 2 + 132£ 3 ) 

+ c 3 f (16 + 114£ + 265£ 2 + 1002£ 3 + 212£ 4 + 776£ 5 ) 

+ c 4 (2 + 16^ + 56£ 2 + 202£ 3 - 147£ 4 + 318£ 5 - 148£ 6 + 116£ 7 ) 



+ c 5 f (-57 + 104£ - 37^ + 68£ d - 24£ 4 ) 
+ c 6 £ (20 - 11£ + 2£ 2 - 14£ 3 ) - c 7 (3 + 2£ + 2£ 2 ) 

+ c (y 5 - y 4 ) (e^ 5 (1 + 8£) + c£ 3 (2 + 33£ + 86£ 2 + 368£ 3 ) 

+ c 2 £ 2 (12 + 82£ + 323£ 2 + 158£ 3 + 536£ 4 ) 

+ c 3 (2 + 16^ + 56 f + 22£ 3 + 313£ 4 - 6£ 5 + 192£ 6 ) 

+ c 4 f (57 - 64£ + 33£ 2 - 52£ 3 + 8£ 4 ) + c 5 £ (-20 + 5£ - 6£ 2 + 10^ 3 ) 



+ 



c 6 (3 + 2£ + 2£ 2 ) ) 



(A.13) 



while Ps(£, c) can be obtained from p^, c) by simply interchanging and 2/5. 

Having obtained J^;(^,c), it is a simple matter to derive the 'hard' component of the 
characteristic function, JF ft (^, c) in Eq. (2.18). In fact, since the only change is in the lower 
limit of integration, which for Th coincides with one of the branch points of the integrand, 
one obtains the same linear combination of elliptic integrals, with each incomplete inte- 
gral replaced by the corresponding complete one. Furthermore, the 'remainder' rational 
function which would play the role of r(£, c) in the present case vanishes. The result is 

Fh&c) = /(£,c)K[m(£,c)] + e (£, c) E \m (£, c) J + p (£, c) II [n (£, c) , m (£, c) 

+ P4 (£ c) n n 4 (£, c) , m (£, c) + p 5 (f, c) n n 5 (f , c) , m (£, c) . (A.14) 

Clearly, the same technology used here can be applied to the simpler situation in which 
£ = 0. The characteristic function in this limit coincides with the leading-order coefficient 
for the C-parameter distribution, which, to our knowledge, was never computed in closed 
form before. In the limit £ — > one finds 



yi 

y2 

y3 



-- + o{f,c 3 ) , 



1 + 4c - VI -8c 

2(1 + c) 
1 + 4c + VI -8c 

2(1 + c) 



+ O(0 , 
+ O(0 , 



(A.15) 



while 7/4^0 and 7/5 — > 1. At this point one can get to the result by either taking the limit 
£ — > in Eq. (2.18), and then applying the algorithm to express that integral in terms of 
the basic set of (complete) elliptic integrals, or one can take directly the limit of Eq. (A.14), 
since in the massless limit only Th contributes. Using either method, the result takes the 
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form 



^b(c) = /o(c) K mo(c) + eo(c) E 



m (cj 



+ p (c)n n (c),m (c) 



(A.16) 



where, as the notation suggests, the various functions are the £ — ► limits of the corre- 
sponding functions in Eq. (A. 14). Some care must be exercised in taking the limit, due 
to the 1/yi singularity in /(£, c). The coefficients £>4(£, c) and ps(£, c) vanish in the mass- 
less limit (while the corresponding integrals are nonsingular) , so that only three complete 
elliptic integrals appear in the final answer. Explicitly one finds 



m (c) 
n (c) 
/o(c) 

eo(c) 
Po(c) 



2^1 -8c 



1 - 4c(l + 2c) + VI - 8c ' 
4V1 - 8c 

(TTTr^W ' 

4V2(l-2c(2 + c)) 



(1 + c) 3 v/ 1 - 4c(l + 2c) + Vl^c 

3(1 + 2c) v/l - 4c(l + 2c) + VT^8c 

V2c(l + c) 3 
V2(2 + c + 2c 2 ) (1 - VT^8^) 2 



;(1 + c) 3 v/l-4c(l + 2c) + VT^ 



8c 



(A.17) 



As shown in Appendix B, Eq. (A.17) reproduces known results for the leading singular 
behavior near c = [13]. 



B. Expansions of the characteristic function 

Here we provide some details on the expansions of the characteristic functions T\ and Th at 
small c. As one can see from Eq. (3.2), and as discussed in Section 3, the relevant scaling 
limits for the computation of log-enhanced contributions to the Sudakov exponent are 
c — ► with a = £/c kept constant for J^i, and c — > with b = £/c 2 kept constant for both 
JFz and Th- These limits can be computed in two different ways: either carefully taking 
the appropriate limit of the integrand in Eq. (2.18), and then performing the resulting 
simplified integration, or by expanding directly the elliptic integrals in Eqs. (A. 8) and 
(A. 14). In either case one needs the expansions of the roots of the cubic equation (2.16). 
For constant a = £/c one finds 



Vi 



V2 



2(1- vTT^)c-2af2-q- 2 " (1 q ) ) c 2 + O (c 3 ) , 

V Vl + a? J 

2(1 + Vl + ^)c -2a f2-a+ 2 ~ / ^Z T Q) ^c 2 + 0(c 3 ) , 

V VI + a 2 / 



y 3 = 1 - (1 - 2a)c - (3 - 6a + 4a 2 )c 2 + O (c 3 ) . (B.l) 
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On the other hand, for constant b = £/c 2 one finds 

Vi = ~b 2 c 3 + O (c 5 ) , 

y 2 = 4c + (4 - b) 2 c 3 + O (c 4 ) , 

y 3 = 1 - c - (3 - 2b)c 2 - (13 - 6b)c 3 + O (c 4 ) . (B.2) 

Similar expansions can immediately be derived for the other singular points, yo, y^ and y§. 

Working at the level of the integrand, f(y, c, £) in Eq. (2.19), the limit c — > with 
a = £/c constant is particularly simple. In that region in fact both limits of integration 
behave like constants, and one can simply expand f(y, c, £) in powers of c at fixed a, 
obtaining 

x 2(2 - (2 + q)j/ + y 2 ) ^ 

f{y,c,ac) = + C(c), (B.3) 

cy 

which can be immediately integrated to give the first line of Eq. (3.3). The limit with fixed 
b = £/c 2 is slightly more difficult, since in that region the lower limits of integration for 
both Ti and Th vanish linearly with c. Simply expanding the integrand in powers of c at 
fixed b is not sufficient in this case, since all powers of c contribute to the leading behavior 
of the integral. The reason is easily tracked to the square-root singularity at y = y 2 . One 
can then solve the problem by keeping exactly the factor (y — y 2 )~ 1 ^ 2 , and expanding the 
other factors of f(y, c, £) in powers of c at fixed b, with the result 



9n 2(2-2y + 2/ 2 ) , n , 
f(y,c,bc 2 ) = V _ +O(c ) 
c^y^y - y2 



(B.4) 



It is fairly easy to see that Eq. (B.4) gives the second line of Eq. (3.3) when integrated 
between be and ys, while it gives Eq. (3.4) when integrated between y 2 and y%. 

The same results were obtained by expanding directly the final expressions for T\ 
and Th in Eq. (A. 8) and Eq. (A. 14), respectively. This involves expanding the elliptic 
functions in the approriate scaling limits, a non-trivial task owing to the singularities 
of these functions. One effective method to perform this expansion is to first scale the 
integration variable such that the integration limits are constants and then express the 
product of square root factors as a single denominator using Feynman parametrization. 
Let us demonstrate this in the case of the incomplete n(n, 4>, m) integral which contributes 
to the C(l/c) terms in T\ in both limits. Since the coefficient multiplying this function in 
c) is 0(c) in these limits, 

p(e,c)^( 32(1 -° )2c fiX6da (B.5) 
; \ 32c fixed b V ; 

relevant terms in II [n(£, c), x(£, c), m(£, c)] would be C(l/c 2 ). Starting from the definition 
(A. 2) and defining z = sincp, one gets 



U(n,<f>,m) = - [ , UW (B.6) 



dw 

y/w (1 — nz 2 w)\J (1 — z 2 w){l — mz 2 w) 



z f 1 da f 1 dw 1 
2vr J y/ay/1 - a J yfw 1-n 



nz 2 w a(l — z 2 w) + (1 — a)(l — mz 2 w) 
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The w integral can be easily done, and the result be safely expanded in the two relevant 
limits. After the expansion the integration over the Feynman parameter a can be readily 
performed. The leading terms are the following: 



n[n(£,c), X (£,c),m(£,c)] 



1 1 



x < 



In a 




fixed a 



fixed b 



(B.7) 



In the case of JF/, in both limits, additional 0(l/c) contributions appear in Eq. (A. 8) from 
£[x(£,c),m(£,c)] =1 + O(c ), since 



e(£,c) 



-(3 + 2a)/c 
-3/c 



fixed a 
fixed 6 



(B.8) 



Finally for fixed a also r (£, c) contributes: r (£, c) ~ a(4+a)/c. Using the method described 
above it is straightforward to verify that the other elliptic integral terms in Eq. (A. 8) do 
not contribute at order 0(l/c). The results are summarized by Eq. (3.3). The case of Th 
is simpler: the elliptic integrals in Eq. (A. 14) are the complete ones and only the fixed b 
limit is relevant. One gets 0(l/c) contributions from U[n(£,c),m(£,c)] ~ — ln(c)/(8c 2 ) 
and from E [m(£, c)] = 1 + O(c ). The result is summarized by Eq. (3.4). 

Finally, we note that using the same techniques one can also treat the exact expression 
for the leading order, which is given by the characteristic function at £ = (see Eq. (A. 16)). 
The relevant elliptic integrals have the asymptotic behavior 



n 



K m (c) 
E m (c) 
n (c),TO (c) 



--logc + O (c 3 logc) , 
l + 0(c 3 logc) , 

- ^ log c - — (1 + log c) + O (c° log c) 



(B.9) 



which leads to 



3 + 4 log- c 

jr Q ( c ) = _ °^ " & +l-281ogc + 0(clogc) 



(B.10) 



in agreement with Ref. [13], when the overall normalization is taken into account. 
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